AD  A 0 3 8 0 5 5 


AFATL-TR-76-122 


STUDIES  ON  THE  PENETRATION 
MECHANICS  OF  EGUN  SAND 


VULNERABILITY  ASSESSMENTS  BRANCH 
WEAPONS  SYSTEMS  ANALYSIS  DIVISION 


.j  D C 

;l  rfFTr  nn/7[? 


OCTOBER  1976 


i W APR  12  1977 

i n i 


uur^gj[f~El 


FINAL  REPORT-JULY  1976  TO  SEPTEMBER  1976 


Approved  for  public  release;  distribution  unlimited 


AIR  FORCE  ARMAMENT  LABORATORY 

AIR  PORCt  SYSTIMS  COMMAND*  UNITID  STATIS  AIR  FORCt 


IGLIN  AIR  FORCE  BASE,  FLORIDA 


This  Document  Contains  Page/s 
Reproduced  From 
Best  Available  Copy 


mmmMM 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  »«St  ('tF'lorl  Dai*  fe'rit mfd) 


0, 


REPORT  DOCUMENTATION  PAGE 


TT  REPORT  JUMPER 


AFATL-TR-76-122 


77 


2,  GOVT  ACCESSION  NO. I 


TITLE  (mnd  SubllJUJ 


STUDIES  ON  THE  PENETRATION  MECHANICS  OF  ) 
2GLIN  SAND  . 7?- ' 


|»  A PFP.r>n  rQYJEaiO^ 


7.  ^AUTHORf»J 

John  77S . lins 
Robert  L.ySierakowski 


/ 

/ 


I.  PERFORMING  ORGANIZATION  NAME  AND  AODRESS 

Vulnerability  Assessments  Branch 
Weapon  Systems  Analysis  Division 
Air  Force  Armament  Laboratory 
Eglin  Air  Force  Base.  Florida  32542 


It.  CONTROLLING  OFFICE  NAME  AND  AODRESS 


Air  Force  Armament  Laboratory 
Armament  Development  and  Test  Center 
Eglin  Air  Force  Base,  Florida  52542 


14.  MONITORING  AGENCY  NAME  ft  ADDRESSf//  Jltteront  from  Controlling  Otllco) 


tC  DISTRIBUTION  STATEMENT  (o I t Mm  Roport) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


J.  PECIPTfT'S  CATALOG  NUMBER 


Final  . Ju 

6 PERFORMING  ORG,  R'fSiRf  nllmbe 


ER 


a.  CONI  RACT  OR  GRANT  NUMBFRfuJ 


10.  PPOGPAM  F.LEMENT,  PROJECT.  TASK* 
APEA  A WORK  UNIT  NUMBERS 


Project  N0/7ST2549 
Task  No. 

Work  Unit  NoV^7  0OT 


'JMBER  OF  PAGES 


29 


15.  SECURITY  CLASS,  (ot  thlm  roport) 

UNCLASSIFIED 


I5«.  DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ol  the  abstract  entered  In  Block  20,  If  different  trom  Report) 


C>  ^ /e  J,,  t~ 


t8.  SUPPLEMENTARY  NOTES 


Available  in  DDC. 


^9.  KEY  WORDS  (Continue  on  revetae  aide  It  necaaeary  and  identity  by  block  number) 


Penetration  Mechanics 
Flash  X-rays 
Poncelet  Prediction 
Sandia  Empirical  Model 


Sand 


20.\aBSTRACT  (Continue  on  revarae  aide  It  necaaeary  and  Identify  by  block  number) 


DD  F0RM 

1 JAN  71 


1473  COITION  OF  1 NOV  65  IS  OBSOLETE 


UNCLASSIFIED 


systematic  study  of  predictive  techniques  for  earth  penetration  has  been 
reviewed.  Based  upon  this  study  an  industrial  proprietary  computer  code,  which 
appears  to  be  the  most  flexible  and  extensible  existing  modal  currently  avail- 
able, has  been  used  to  examine  controlled  experimental  data  generated  at  AFATL. 
Results  of  this  data  analysis  indicate  that  use  of  the  code  requires  some  judg- 
ment in  proper  selection  of  input  parameters.  In  addition,  test  data  generated 
at  AFATL  indicate  that  the  drag  coefficient  in  predictive  soil  penetration  


SECURITY  CLASSIFICATION  OF  THIS  PAGE  flWi»n  DtIM  Enlorod) 


U (TO 


■ran 


T 


SECURITY  CLASSIFICATION  OF  THIS  P *GE(Whvn  Dura  Enfrvd) 


PREFACE 


This  report  documents  work  performed  during  the  period  from  July  1976 
to  September  1976  by  the  Vulnerability  Assessments  Branch,  Weapon  Systems 
Analysis  Division,  Air  Force  Armament  Laboratory,  Eglin  Air  Force  Base, 
Florida  32S42,  under  Project  2549.  Messrs.  John  A.  Collins  and  Robert  L. 
Sierakowski  (DLYV)  conducted  the  study  for  the  Armament  Laboratory. 

This  report  has  been  reviewed  by  the  Information  Office  (01)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS) . At  NTIS, 
it  will  be  available  to  the  general  public,  including  foreign  nations. 


This  report  has  been  reviewed  and  is  approved  for  publication. 
FOR  THE  COMMANDER 

(A  A ytjuAAM  W 

J./R.  MURRAY  // 

Chief,  Weapon  Systems  /^yflysis  Division 


i 

(The  reverse  of  this  page  is  blank) 


TABLE  OF  CONTENTS 


Section  Title  Page 

I INTRODUCTION  1 

Purpose  and  Scope  1 

Background 1 

II  EXPERIMENTAL  DETAILS 3 

Test  Procedure 3 

III  THEORETICAL  CONSIDERATIONS  9 

Description  of  Model  9 

Penetration  Predictions  . 11 

IV  DISCUSSION 17 

V CONCLUSIONS  20 

References  21 


iii 


T.  1ST  OF  FIGURES 


Figure  Title  Page 

1 Blunt-Ended  Cylindrical  Projectiles  4 

2 Experimental  Test  Set-Up  5 

3 Scaled  X-Ray  of  Shot  Number  26  6 

4 Schematic  of  Coil  Positions  8 

5 Coordinate  Frames  of  Reference  . . 10 

6 Schematic  of  X-Ray  Positions  12 

LIST  OF  TABLES 

Table  Title  Page 

1 Experimental  Test  Matrix  7 

2 Comparison  of  Experimental  Data  and  Predictions  - 

Dry  Sand 14 

3 Comparison  of  Experimental  Data  and  Predictions  - 

Wet  Sand 15 

4 Drag  Coefficients  Versus  Velocity  - Wet  Sand  ....  19 

5 Drag  Coefficients  Versus  Velocity  - Dry  Sand  ....  19 


LIST  OF  ABBREVIATIONS,  ACRONYMS  AND  SYMBOLS 


Elemental  Force  (kg) 

Elemental  Area  (m2) 

Medium  Density  (kg/m^) 

Wave  Speed  of  Target  Medium  (m/sec) 

Local  Absolute  Velocity  Vector  (m/sec) 

Exponential  Decay  Factor  (1/sec) 

Time  (sec) 

Dummy  Variable  Time  (sec) 

Normal  Flow  Coefficient 

Shear  Flow  Coefficient 

A Function  of  Medium  Density  and  Local  Incidence  Angle 
Medium  Resistance  Coefficient  (kg/m2) 


v 

(The  reverse  of  this  page  is  blank) 


SECTION  I 


INTRODUCTION 


PURPOSE  AND  SCOPE 

The  purpose  of  this  study  was  to  use  a limited  distribution  code  (spe- 
cifically, an  industrial  proprietary  code  developed  by  AVCO  Corporation)  to 
examine  the  extent  of  its  predictive  capabilities  for  comparison  to  exper- 
imental data  generated  at  AFATL.  In  addition,  because  of  the  apparent  com- 
prehensiveness of  the  code  (three  dimensional,  six-degrees-of-freedom  model) 
in  characterizing  penetrator  and  target,  a systematic  screening  of  penetration 
variables  was  undertaken  to  obtain  quantitative  information  on  the  importance 
of  various  properties  on  penetrator  performance. 

BACKGROUND 

Early  studies  on  determining  the  penetration  mechanics  of  impacted  media 
have  been  limited  to  predictions  of  the  depth  of  penetrators  into  target 
media,  in  recent  years  more  advanced  analytical  models  have  evolved  which 
predict  target/penetrator  interaction  as  well  as  the  penetrator  trajectory. 
Although  the  models  introduced  appear  to  have  considerable  merit,  verifi- 
cation of  their  usefulness  has  remained  somewhat  unchallenged  due  to  the 
limited  experimental  data  available  for  predictive  correlation. 

In  general,  the  theories  advanced  have  been  categorized  in  a recent  survey 
(Reference  1)  as  being  of  the  semianalytical , analytical,  theoretical,  and 
empirical  type.  The  semianalytical  models  include  the  earliest  penetration 
equations  proposed  (References  2,  3,  4)  and  rely  to  a large  extent  on  exper- 
imental data  for  evaluation  of  equation  constants.  The  analytical  techniques 
proposed  are  extensions  of  the  one-dimensional  semianalytical  models  and 
introduce  added  target/penetrator  constitutive  equations  for  better  char- 
acterizing penetrator  performance.  Principal  among  these  models  are  the 
cavity  expansion  theory  (Reference  5)  and  the  differential  area  force  law 
(Reference  6) . The  cavity  expansion  theory  characterizes  the  target  material 
compressibility  using  an  elastic-plastic  working  material  model.  This  ana- 
lytical technique  has  been  extended  (Reference  7)  to  include  such  variables 
as  penetrator  nose  shape,  target  layering,  and  impacts  deviating  at  small 
angles  relative  to  the  target  surface  normal.  The  range  of  applicability 
of  such  models,  however,  appears  to  be  limited  to  essentially  normal  impact 
configurations.  An  apparently  more  versatile  analysis  which  is  based  upon 
a three-dimensional,  six-degrees-of  freedom  model  and  which  includes  oblique 
impact  has  been  introduced  in  Reference  6.  This  analytical  model  considers 
the  penetrator  as  a rigid  body  for  establishing  the  loads  and  accounts  for 
pressure  induced  drag,  target  structural  resistance,  and  penetrator  surface 
effects.  Once  the  loads  are  established  by  this  code,  they  can  be  used  in 
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model  analysis  models  or  equivalent  techniques  to  estimate  loads  transmitted 
to  interior  penetrator  components.  Theoretical  techniques  (References  8,  9, 
10)  rely  upon  formulating  constitutive  equations  of  the  target  and  penetrator 
materials  for  prediction  of  the  governing  equation  parameters.  Many  of  these 
techniques,  as  currently  available,  are  extensions  of  target  penetration  codes 
based  upon  finite  difference  or  finite  element  methods,  and  generally  consider 
symmetric  vehicle  shapes.  While  allowing  for  continuous  deformation  of  both 
target  and  penetrator,  such  computer  models  involve  accurate  modelling  of  the 
penetrator  and  target  properties  for  characterizing  penetrator  performance. 
Finally,  empirical  predictions  of  penetrator  performance  such  as  those  in 
Reference  11  have  found  usefulness  for  specific  velocity  regimes  and  vehicle 
sizes.  Close  examinatiqn  of  such  predictor  equations  shows  that  many  of 
these  formulations  are  modified  forms  of  the  one-dimensional  classical  pene- 
tration equations  (References  2,  3,  4). 
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SECTION  II 


EXPERIMENTAL  DETAILS 


TEST  PROCEDURE 

The  experimental  projectile  penetration  data  recorded  has  been  obtained 
at  the  AFATL  test  range  using  a facility  specifically  converted  for  pene- 
tration type  testing.  The  test  facility  has  provisions  for  adaptation  so 
that  both  horizontal  and  vertical  firings  can  be  accommodated. 

In  the  present  experiments,  a series  of  0.02  by  0.23  meter  long,  blunt- 
ended  cylindrical  steel  projectiles  (AISI-W1)  (Figure  1)  were  fired  hori- 
zontally into  a 1.2-meter  long  by  0.15  meter  by  0.40  meter  open  top  box 
(Figure  2).  The  test  chamber  was  slowly  filled  with  Eglin  sand  sieved  to 
remove  large  debris,  using  a U.S.  Standard  Sieve  Series  Number  25  screen. 

A 20 millimeter  gun  was  used  for  firing  the  projectiles  with  velocity  con- 
trolled by  the  powder  charge  used.  For  the  present  tests,  these  velocity 
regimes  have  been  examined  for  flat-nosed  projectiles  fired  into  dry  and 
wet  sand  targets.  The  initial  impact  angle  of  incidence  was  approximately 
zero  degrees  for  all  tests.  A summary  of  the  specimens  tested,  with  target 
conditions  indicated,  is  included  in  Table  1. 

In  order  to  obtain  data  for  examining  penetrator  trajectory  through  the 
sand  medium,  vehicle  stability,  and  forces  acting  on  the  projectile,  a 
number  of  sensing  elements  were  tried  and  from  these  several  were  selected 
for  recording  the  test  data.  The  most  successful  of  the  data  collection 
methods  used  was  found  to  be  flash  radiography.  For  the  current  experiments 
five  sequentially  spaced  X-ray  heads  were  used.  The  first  150  KV  unit  was 
located  0.038  meter  from  the  front  of  the  box  while  the  remaining  four  300  KV 
units  were  spaced  0.38  meter  on  center  from  one  another.  The  units  were 
positioned  at  nominal  standoff  distances  of  0.55  meter  from  the  box  center 
line.  To  determine  the  position  of  the  vehicle  down  the  test  chamber,  a 
series  of  letters  (A  through  Q)  were  attached  in  a horizontal  line  along  the 
wall  of  the  box.  The  letters  were  spaced  at  nominal  intervals  of  0.07  meter 
with  respect  to  one  another  and  at  a distance  of  0.20  meter  from  the  top  of 
the  box.  These  letters  provided  a ready  means  of  locating  the  horizontal 
and  vertical  position  of  the  penetrator  in  the  X-ray  record  for  any  instant 
of  time  during  penetrator  motion  through  the  box.  In  order  to  detect  if  any 
changes  in  soil  motion  had  occurred  during  passage  of  the  projectile  through 
the  box,  a series  of  0.0015  meter  steel  markers  were  suspended  in  the  soil 
medium.  Both  preshot  and  post  shot  X-ray  records  were  made  to  provide  this 
information  with  Dupont  Lightning  Plus  being  used  as  the  X-ray  film.  An 
overview  of  the  test  chamber,  X-rays,  and  film  attached  in  a test  setup  is 
shown  in  Figure  2,  while  a scaled  schematic  of  the  X-ray's  positions  in  a 
horizontal  plane  is  shown  in  Figure  3. 
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TABLE  1.  EXPERIMENTAL  TEST  MATRIX 


Velocity 

Projectile 

Type 

Target 

210  M/Sec 

320  M/Sec 

400  M/Sec 

Flat  Nosed 

Dry  Sand 

15  16  17  18  19 

20  22  23  24 

14  25  26  27  29 

Flat  Nosed 

Wet  Sand 

70  71  72  73 

36  37  38  74  75 

76  81  82  83  84 

As  mentioned,  projectile  striking  velocites  were  controlled  by  varying 
the  powder  load  in  a primed  20mil limeter  case  with  striking  velocity  measured, 
using  paper-back  velocity  screens  located  at  fixed  intervals  near  the  test 
chamber  entrance.  The  X-ray  units  were  triggered  by  the  use  of  a foil  make 
switch  with  timing  sequence  predicted  according  to  best  estimates  furnished 
by  the  Project  Engineer. 

In  addition  to  the  X-ray  units,  a series  of  round  copper  coils  were 
attached  to  the  walls  of  the  test  chamber  at  fixed  intervals  for  use  as  a 
further  check  on  projectile  velocity  while  in  the  test  chamber.  Such  measure- 
ments were  made  possible  by  magnetizing  the  projectiles  before  each  test  with 
a nominal  field  strength  of  150  Gauss  recorded  at  the  front  and  tail  section 
of  each  penetrator.  Static  measurements  were  made  of  projectile  motion 
through  the  center  of  the  coil  to  establish  peak  magnetic  response  of  the 
projectile  while  cutting  the  flux  lines.  Deviations  in  peak  magnetic  re- 
sponse were  measured  in  the  order  of  0.02  to  0.04  meter  in  from  the  nose 
and  tail  of  the  projectile  and  this  range  was  considered  the  error  bound 
for  the  dynamic  tests.  The  size  of  the  coils  used  (0.15  meter  diameter) 
was  determined  by  trying  to  insure  that  the  magnetic  lines  of  force  would 
be  cut  at  an  optimum  position  which  presumed  that  the  penetrator  path  was 
essentially  straight  and  true  through  the  center  of  the  test  chamber.  For 
the  current  vehicle  shapes  and  velocity  regimes  this  was  found  to  be  essen- 
tially true  with  errors  due  to  deviations  in  straight  line  trajectory  and 
variable  location  of  peak  magnetism  on  the  projectile  to  be  of  the  order  of 
five  percent.  A schematic  of  the  coil  positions  used  for  the  tests  reported 
herein  is  shown  in  Figure  4 . 


SECTION  III 


THEORETICAL  CONSIDERATIONS 


DESCRIPTION  OF  MODEL 

Of  the  current  methods  available,  the  AVCO  developed  code  known  in  the 
literature  as  the  differential  area  force  law  (DAFL)  appears  to  be  the  most 
versatile  and  represents  the  only  one  which  can  handle  non-normal  impacts 
and  .distributed  longitudinal  as  well  as  circumferential  surface  loadings. 
Although  not  well  documented  in  the  literature,  it  appears  to  merit  con- 
sideration based  upon  its  apparent  predictive  capabilities. 

The  form  of  the  DAFL  model  described  here  constitutes  an  extension  of 
previous  one-  and  two-dimensional  developments  and  includes  the  complete 
three-dimensional  simulation  of  projectile  penetration.  The  three  trans- 
lational and  three  rotational  equations,  which  describe  the  penetrator 
motion,  rely  upon  knowledge  obtained  from  experiment  of  certain  penetration 
properties  and  consider  normal  as  well  as  non-normal  impact  simulation.  As 
previously  mentioned,  such  types  of  models  have  been  classified  as  semi- 
analytical.  The  essential  features  of  this  model  are  outlined  in  References 
6,  12,  13,  14,  15  and  16  and  a brief  description  is  included  below. 

Four  coordinate  frames  of  reference  are  used  in  the  analytical  model. 
These  are  a fixed  inertial  frame  in  the  target  medium,  a target  reference 
frame,  a body  fixed  reference  frame  at  the  center  of  gravity  of  the 
projectile,  and  a pseudo-body  reference  frame  passing  through  the  center 
of  gravity  of  the  projectile.  For  the  initial  computations,  the  projectile 
is  considered  as  a rigid  body  moving  in  a fluid  medium  (sand)  which  can  be 
subdivided  into  a series  of  surface  elements  along  the  circumference  of 
the  disk  section  (Figure  5) . The  total  number  of  disk  sections  considered 
at  one  time  is  less  than  200,  while  the  number  of  sector-like  segments  is 
greater  than  six.  The  sections  and  sectors  together  delineate  an  elemental 
surface  area  on  which  the  elemental  force  acts.  This  elemental  force  is 
considered  to  be  made  up  of  the  following  components:  (1)  Normal  and 

shear  flow  force  terms  proportional  to  velocity  squared  and  defined  in 
terms  of  form  and  surface  drag  coefficients,  (2)  a transient  shock  term 
which  depends  upon  the  acoustic  impedance  of  the  target  medium,  (3)  a 
cratering  term  which  accounts  for  the  unequal  pressure  distribution  along 
the  length  of  the  projectile,  (4)  a chipping  term  which  accounts  for  unequal 
pressure  distribution  around  the  circumference  of  the  penetrator,  and  (5) 
a friction  term  related  to  the  static  bearing  strength  resistance  of  the 
target  medium.  Expressed  analytically  the  DAFL  is  written  as: 

f "at  2 

dF  = dA  [_  pcve  u(t  - T)  + Cnkv 

+ CT  kv2  + n J + Target  Surface  Effects 


9 


Each  of  the  elemental  force  components  acting  on  the  respective  area 
elements  is  transferred  to  the  center  of  gravity  of  the  vehicle  and  summed 
to  produce  a total  resultant  force  along  the  body  axis.  Moments  of  the 
differential  forces  are  calculated  about  the  center  of  gravity  of  the  pro- 
jectile without  the  necessity  of  calculating  a center  of  pressure.  Thus, 
the  necessity  for  introducing  other  aerodynamic  coefficients  in  addition  to 
the  normal  and  tangential  drag  coefficients  incorporated  in  the  model.  Such 
considerations  would  appear  to  enhance  the  physical  completeness  of  the  model 
and  its  corresponding  predictive  capabilities.  Specifically,  less  ambiguity 
would  be  introduced  into  assessing  the  manner  of  input  accociated  with  the 
cratering  and  chipping  terms. 

The  necessity  for  experimental  data  as  input  to  the  force  law  becomes 
readily  apparent  in  the  parameters  of  the  terms  appearing  in  the  analytical 
expression.  (See  List  of  Abbreviations,  Acronyms  and  Symbols  for  term  iden- 
tification.) The  computer  code  handles  certain  standard  type  nose  shapes 
through  attached  subroutines;  however,  other  shape  types  can  be  handled  by 
substituting  specific  input  information. 

PENETRATION  PREDICTIONS 

As  discussed  in  the  experimental  section,  data  on  the  velocity  trajectory 
of  20millimeter  projectiles  in  dry  and  wet  sand  have  been  obtained  at  AFATL. 
With  the  use  of  X-ray  radiography  the  trajectory  of  the  impactor  can  be  mon- 
itored sequentially  with  corresponding  angle  of  attack  of  the  projectile  re- 
corded. A typical  example  of  the  data  collected  by  this  technique  is  shown 
in  Figure  6 which  documents  Test  Number  26  from  Table  1.  Thus,  the  data  as 
obtained  from  these  experiments  provides  and  independent  basis  (other  than 
AVCO’s)  for  testing  the  DAFL  code  for  predicting  the  experimental  results. 

In  addition,  the  experimental  data  obtained  coupled  with  the  code  predictive 
capability  provides  a rational  basis  for  evaluating  the  significance  of 
individual  contributions  to  the  total  resultant  force. 

To  initiate  data  reduction  between  stations  for  each  shot,  the  tests 
indicated  in  Table  1 were  catalogued  and  data  on  the  striking  velocity’, 
projectile  position,  and  time  of  arrival  at  the  corresponding  position 
recorded.  This  data  was  corrected  for  any  photographic  distortion  in  the 
plane  of  the  projectile  by  scaling  the  true  vehicle  dimensions  with  respect 
to  the  photographically  recorded  size.  This  adjusted  data  was  then  used  to 
reevaluate  and  record  an  adjusted  projectile  position  with  subsequent  anal- 
ysis of  vehicle  velocity,  attitude,  and  other  parameters  evaluated.  Some 
additional  comments  related  to  this  data  reduction  procedure  are  included 
in  Reference  17.  This  report  also  includes  a comparison  between  tests  of 
the  one-dimensional  form  drag  flow  term  in  the  Poncelet  equation  versus  the 
corresponding  empirical  model  of  Sandia. 

The  aforementioned  data  was  then  used  in  a Poncelet  predictive  equation 
to  obtain  information  on  the  drag  coefficient  as  a function  of  velocity 
which  in  turn  could  then  be  used  to  determine  an  arbitrary  displacement  versus 


velocity  and  time  profile  for  each  shot.  A summary  of  this  data  is  included 
in  Tables  2 and  3 for  all  of  the  test  cases  indicated  in  Table  1.  Since  a 
spectrum  of  data  is  indicated  in  Tables  2 and  3 for  each  shot,  and  particulars 
within  each  shot,  it  was  possible  to  establish  a drag  coefficient,  versus 
velocity  profile.  This  information  has  importance  for  comparison  with  avail- 
able DAFL  input  information  in  Cp,  C versus  velocity  which  is  described  below. 

In  order  to  evaluate  the  codes  predictive  capabilities,  an  examination 
of  one  of  the  AFATL  tests  described  in  Table  1 (Number  26)  was  used  as  a 
first  reference.  To  avoid  introduction  of  any  input  bias,  the  code  was  run 
using  the  standard  table  of  lookup  functions  and  tabulated  data  for  sand  as 
described  in  Reference  12.  Initial  results  obtained  from  the  code  did  not 
provide  satisfactory  agreement  with  the  experimental  data.  Subsequently, 
an  examination  of  input  data  for  each  term  was  used  to  systematically  screen 
the  influence  on  and  contribution  to  the  total  resultant  force.  As  a con- 
sequence of  the  results  obtained  from  the  code  and  for  the  velocity  regime 
selected  for  study,  it  was  determined  that  the  flow  terms  constituted  the 
dominant  terms  in  the  force,  particularly  the  drag  term.  This  result  had 
led  to  a reexamination  of  the  DAFL  code  input  data  and  to  locating  an 
apparent  inconsistency  in  the  table  of  lookup  function  data  for  C^,  CT 
versus  velocity. 

On  the  basis  of  the  AFATL  test  data  a Cp  versus  velocity  lookup  function 

table  was  constructed  and  inputed  into  the  DAFL  program.  In  addition,  the 
DAFL  code  was  run  using  then  established  table  of  lookup  function  data 
supplied  with  the  code  from  Picatinny  Arsenal  for  Cp  versus  velocity  with 

all  other  input  parameters  held  constant  between  the  two  runs.  A summary  of 
the  results  obtained  from  these  runs  is  included  in  Tables  2 and  3.  As  can 
be  seen,  the  DAFL  program  run  with  Cp  data  extrapolated  from  the  AFATL 

program  provides  for  very  good  agreement  between  code  and  experimental  data. 

On  the  other  hand,  the  AVCO  tabulated  data  is  shown  to  provide  rather  un- 
realistic predictions  on  both  depth  of  penetration  and  resultant  velocity. 
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TABLE  2.  COMPARISON  OF  EXPERIMENTAL  DATA  AND  PREDICTIONS  - DRY  SAND 


Shot 

No. 

Striking 

Velocity 

(m/sec) 

(•?) 

Drag 

(cd) 

Distance 

Cm) 

Avg  Vel 
(m/sec) 

Cal  Vel 
(m/sec) 
AVC0 

Cal  Vel 
(m/sec) 
AFATL 

at  0.16  «= 

16 

212.4 

0.7621 

1.756 

0.2723 

178.7 

94.9 

171.3 

0.8495 

1.958 

0.5154 

141.3 

140.3 

17 

212.1 

0.772 

1 . 779 

0.3198 

174.1 

94.8 

164.6 

0.690 

1.590 

0.5898 

139.3 

131.7 

0.889 

2,0^9 

0.9009 

108.6 

101.2 

18 

213.4 

0.819 

1.887 

0.3068 

174.6 

95.0 

167.3 

0.726 

1.673 

0.5700 

139.3 

134.7 

1.047 

2.413 

0.8260 

106.5 

108.6 

19 

210.6 

0.802 

1.849 

0.3109 

176.5 

164.5 

0.816 

1.880 

0.5728 

139.6 

132.5 

0.932 

1.896 

0.8799 

106.5 

102.2 

at  0.21  = 

20 

329.2 

0.728 

1.678 

0.2949 

268.1 

96.0 

263.2 

0.829 

1.910 

0.5519 

214.8 

216.4 

0.732 

1.687 

0.8621 

168.3 

169.4 

23 

328.2 

0.874 

2.014 

0.2990 

269.5 

261.6 

0.684 

1.576 

0.5621 

216.0 

214.1 

0.934 

2.152 

0.8631 

168.7 

168.7 

24 

327.4 

0.713 

1.643 

0.2690 

265.9 

266.9 

0.770 

1.774 

0.5250 

217.8 

219.7 

0.784 

1.814 

0.8260 

172.0 

173.4 

25 

406.0 

0.905 

2.086 

0.2769 

329.1 

328.3 

0.688 

1.400 

0.5339 

264.2 

269.8 

0.794 

1.615 

0.8321 

211.3 

214.4 

at  0.23  = 

26 

406.3 

0.774 

1.783 

0.2799 

329.6 

94.1 

328.5 

0.764 

1.760 

0.5380 

266.5 

270.1 

0.657 

1.515 

0.8400 

214.0 

214.6 

27 

408.7 

0.829 

1.909 

0.2809 

331.4 

330.3 

0.749 

1.725 

0.5380 

266.2 

271.7 

29 

40S.0 

0.839 

1,933 

0.2809 

330.8 

327.3 

0.764 

1.761 

0.5220 

266.3 

274.0 

0.705 

1.625 

0.8390 

214.0 

214.1 

TABLE  3.  COMPARISON  OF  EXPERIMENTAL  DATA  AND  PREDICTIONS  - WET  SAND 


Shot 

No. 

Striking 

Velocity 

(m/sec) 

B1 

r:°) 

Distance 

(m) 

Avg  Vel 
(m/sec) 

Cal  Vel 
(m/sec) 
AVC0 

Cal  Vel 
(m/sec) 
AFATL 

70 

209.1 

1,018 

1.7593 

0.2741 

171.7 

171.5 

0.862 

1.4897 

0.4971 

133.6 

141.7 

2.0412 

3.5276 

0.7590 

94.4 

108.7 

71 

207.9 

1.254 

2.1671 

0.2680 

182.3 

171.1 

0.903 

1.5605 

0.5189 

134.2 

137.7 

1.112 

1.921 

0.8260 

98.8 

99.8 

72 

214.0 

0.7539 

1.302 

0.2850 

181.1 

174.9 

0.9035 

1.5614 

0.5410 

142.3 

140.7 

0.9681 

1.6730 

0.8649 

104.9 

100.5 

73 

212.8 

0.816 

1.4102 

0.2870 

180.9 

173.4 

0.745 

1.2875 

0.5489 

143.6 

138.5 

0.9818 

1.6967 

0.8781 

108.3 

98.0 

36 

326.2 

0.685 

1.1838 

0.3211 

279.0 

285.4 

0.517 

0.893 

0.6060 

231.2 

248.2 

37 

336.5 

0.656 

1.1336 

0.2979 

279.6 

298.1 

0.477 

0.824 

0.5591 

237.0 

264.5 

0.705 

1.2183 

0.8829 

196.6 

221.7 

38 

333.1 

0.597 

1.0317 

0.3000 

284.0 

294.8 

0.535 

0.9245 

0.5870 

240.2 

259.6 

0.577 

0.997 

0.9169 

200.1 

215.3 

74 

334.0 

0.718 

1.2403 

0.3091 

282.7 

294.4 

0.534 

0.9228 

0.5799 

233.4 

258.5 

0.712 

1.2304 

0.9169 

189.7 

214.6 

76 

406.0 

0.829 

1.4326 

0.3091 

348.2 

358.4 

0.379 

0.6549 

0.5799 

291.2 

321.3 

81 

333.8 

0.592 

1.0229 

0.3091 

278.7 

294.1 

0.5121 

0.8850 

0.5700 

236.6 

260.5 

0.757 

1.3082 

0.8691 

193.6 

217.2 

82 

404.8 

0.571 

0.986 

0.3081 

343.2 

357.5 

0.373 

0.6446 

0.5621 

299.3 

322.6 

0.289 

0.499 

0.8890 

269.2 

281.2 
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TABLE  3.  COMPARISON  OF  EXPERIMENTAL  DATA  AND  PREDICTIONS 
WET  SAND  (CONCLUDED) 


mm 


Shot 

No. 

Striking 

Velocity 

(m/sec) 

B1 

Drag 

(cd) 

Distance 

(m) 

Avg  Vel 
(m/sec) 

Cal  Vel 
(m/sec) 
AVCO 

Cal  Vel 
(m/sec) 
AFATL 

83 

419.4 

0.711 

1.228 

0.3071 

347.0 

370.6 

0.3516 

0.6972 

0.S471 

298.5 

336.4 

0.664 

1.1475 

0.8539 

258.7 

297.0 

84 

405.7 

0.566 

0.978 

0.3081 

333.3 

358.3 

0.338 

0.S841 

0.5441 

294.0 

325.7 

0.765 

1.322 

0.8430 

253.0 

287.8 

m. 


SECTION  IV 


DISCUSSION 


This  study  was  undertaken  to  enhance  the  Air  Force  in-house  capability 
to  realistically  predict  the  terradynamic  trajectories  of  projectile  pene- 
tration into  soil/concrete.  Toward  this  end,  development  of  a test  facility 
from  which  detailed  experimental  data  could  be  gathered  was  assimilated. 

The  principal  experimental  tool  used  for  obtaining  this  information  was 
flash  radiography.  The  DAFL  three-dimensional  code  as  supplied  by  Picatinny 
Arsenal  was  then  run  in  conjunction  with  established  one-dimensional  pre- 
dictors to  simulate  the  projectile  trajectories. 

The  use  and  running  of  the  three-dimensional  DAFL  code  at  AFATL  has 
resulted  in  numerous  problems.  First,  the  code  had  been  received  from 
Picatinny  Arsenal  in  an  interim  format,  that  is,  before  Picatinny  Arsenal 
believed  the  code  ready  for  general  release.  As  such,  the  code  was  in  a 
nondocumented  form  without  the  enhancement  of  any  comment  cards  and  physical 
explanation.  The  limited  user  reports,  as  available  on  the  DAFL  code,  also 
failed  to  provide  adequate  information  or  proper  checks  for  running  the  code. 
Nevertheless,  with  the  path  trajectory  data  as  reduced  from  the  experimentally 
obtained  X-ray  data,  the  AFATL  was  in  a position  to  independently  examine 
the  DAFL  code  as  well  as  other  predictors. 

As  can  be  seen  from  Tables  2 and  3,  the  DAFL  code  predictions  did  not 
match  the  time-position-velocity  history  data  from  experiments.  It  was 
suggested  that  this  discrepancy  might  be  due  to  inexperience  in  running  the 
code  which  has  been  indicated  to  require  considerable  engineering  judgment. 
This  could  indeed  be  the  case  as  there  is  a lack  of  documentation  and  ex- 
planation of  the  code  as  pointed  out,  for  example,  in  References  1 and  16. 

For  this  reason,  a careful  study  of  the  DAFL  code  was  initiated.  Each  com- 
ponent which  went  into  the  total  force  term  for  each  of  the  modelled  pene- 
trator's  elemental  areas  was  catalogued.  In  checking  the  influence  of‘  the 
various  components  of  the  total  force  for  the  experimental  data  available, 
each  term  of  the  total  force  was  systematically  compared  by  suppressing 
individual  component  terms.  For  the  case  of  sand,  the  shock  term,  the 
medium  resistance  pressure,  the  shear  drag  and  target  surface/edge  effects 
were  found  to  be  small  in  comparison  to  the  fluid  dynamic  pressure  effect 
denoted  in  terms  of  CN»  This  latter  term  was  found  to  be  the  principal 

force  component  responsible  for  position  dependent  deceleration  in  the  sand 
targets  tested.  This  appears  to  be  logical  as  the  projectile  was  fired 
into  the  sand  horizontally  so  that  the  influence  of  soil  stratification, 
soil  separation  and  reattachment,  and  the  other  principal  DAFL  terms  would 
not  be  of  concern  at  these  test  velocities. 

Using  the  input  data  for  CN  as  given  in  the  code,  it  was  found  that  the 
DAFL  code  predicted  the  projectile  velocity  to  slow  down  abruptly  and  then 


after  it  reached  approximately  one-fifth  of  its  initial  value  to  slowly  tend 

to  zero.  AVCO  has  stated  that  values  of  can  vary  from  2.0  at  zero  velocity 

to  4.0  at  very  high  velocities  (Reference  15).  In  addition,  AVCO  states  that 
for  purely  elastic  impact  of  profile  penetrator  interaction,  where  the  re- 
bound velocity  is  equal  to  the  impact  velocity,  takes  on  the  value  4.0 

while  C^,  — 0 . A second  limiting  case  is  said  to  occur  when  = 2.0,  C - 2.0 

which  represents  simply  purely  inelastic  momentum  transfer.  The  results  of 

using  this  data  are  shown  in  Tables  2 and  3 under  the  seventh  column.  A 
lack  of  data  in  the  seventh  column  indicates  that  when  the  code  stopped 
within  the  computer  time  limits  set  forth  the  vehicle's  position  (depth)  was 
less  than  that  at  which  the  first  X-ray  data  was  obtained. 

On  the  other  hand,  the  AFATL  test  data  also  shown  in  Tables  2 and  3 
(second,  third  and  fourth  columns)  has  CD/C  varying  between  2 at  high 
velocities  (inelastic  momentum  transfer) uto  approximately  4 at  zero  velocity 
(elastic  momentum  transfer)  for  dry  sand.  The  inelastic  momentum  transfer 
occurs  at  the  higher  velocities  and  is  accompanied  by  the  formation  of  a 
soil  nose  when  the  velocity  is  high  enough.  The  values  of  CD  (Tables  4 and 

5),  as  determined  by  AFATL,  have  also  been  found  to  agree  with  those  reported 

in  References  18  and  19.  When  these  values  of  CL  (C..)  were  used  in  the  DAFL 

D N 

code,  that  is,  varying  CD  from  2.0  at  high  velocities  to  4.0  at  low  velocities 

in  stepwise  increments,  very  good  agreement  was  found  between  experiment  and 
DAFL  code  predictions. 

Reference  15  includes  a discussion  on  an  updated  version  of  the  force 
law  coefficients  and  assigns  a new  variable  to  define  an  equivalent  fluid 

flow  coefficient,  specifically, 

t 2 2 

c’  = Cv,  sin  6 + C cos  e 
D N t 

where  0 is  the  average  nose  cone  half  angle.  In  the  present  AFATL  version 
of  the  DAFL  code,  the  influence  of  0 on  the  corresponding  and  CT  as  a 

function  of  velocity  must  first  be  calculated  before  entering  this  data  into 
the  code  in  the  table  of  lookup  functions  for  CN  (v)  and  CT  (v) . When  the 

velocity  vector  and  the  axis  of  the  projectile  are  non-coincident,  the  in- 
fluence of  the  corresponding  nose  shape  factor  may  be  of  some  consequence. 

For  the  present  AFATL  experimental  data  considered  in  this  report,  CN  and  CT 

have  been  calculated  for  a specific  nose  shape,  with  the  corresponding  input 
data  related  to  that  case.  Specifically,  0 - 90“  and  therefore  is  equal 
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TABLE  4.  DRAG  COEFFICIENTS  VERSUS  VELOCITY  - WET  SAND 


V = 

CD 

AFATL 

0 

6 

12.7 

4 

25.4 

3 

76.2 

2.2 

152.4 

1.5 

228.6 

1.0 

304.8 

0.7 

2540. 

0.6 

• 

TABLE  5. 

DRAG  COEFFICIENTS 

VERSUS  VELOCITY  - 

DRY  SAND 

V - 

0 

12.7 

25.4 

76.2 

152.4 

228.6 

304.8 

2540. 

CD 

4 

3 

2.2 

2.0 

1.9 

1.75 

1.75 

1.7 

AFATL 

The 

shock  term  analytically  represented  by 

pcVu  (t  - 

x)  exp  (-at)  was 

not  used  in  this  form  in  the  DAFL  code  as  supplied  to  the  AFATL.  The  specific 
form  of  the  shock  term  used  in  DAFL  is  kv  [1  + f£]  u (t  ■ t)  exp  (-a)  where 

the  independent  variable  t has  been  replaced  by  the  velocity  v allowing  a to 
be  considered  constant  with  respect  to  impact  velocity.  Indeed,  as  in  the 
preceding  case  of  the  nose  factor  correction,  the  shock  term  must  be  first 
evaluated  to  obtain  the  input  data  for  shock  as  a function  of  velocity. 

When  this  is  done,  the  input  data  becomes  a special  case  for  Vo  or  the 
striking  velocity. 

Further,  it  is  observed  that  the  DAFL  shock  term  denoted  by  FETA  1 is 
handled  within  the  code  as  a tabulated  lookup  function  as  opposed  to  a 
function  of  the  projectile  striking  velocity.  This  causes  problems  in  in- 
cluding the  influence  of  a shock  term  in  the  DAFL  code  runs  as  related  to 
computing  the  velocity  as  a function  of  the  depth  of  penetration.  From 
Tables  2 and  3,  it  can  be  seen  that  the  experimentally  computed  Cp  generally 

starts  out  high,  decreases,  and  then  increases.  It  would  be  advantageous 
to  describe  the  initial  penetration  by  both  a drag  term  and  a shock  term. 

The  form  of  this  shock  term  should  be  of  a general  form  for  any  target  ma- 
terial and  should  also  incorporate  functional  dependence  on  velocity  and 
exhibit  a characteristic  penetration  time.  Such  modifications  to  the  DAFL 
code  as  mentioned  above  will  be  undertaken  ir.  a future  study. 

Another  area  which  will  receive  additional  study  in  the  future  is  the 
low  velocity  regime.  In  this  regime,  the  n term  or  target  medium  resistance 
and  the  medium  resistance  friction  as  the  soil  reattaches  to  the  projectile 
remain  the  major  influences  on  vehicle  deceleration. 
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SECTION  V 


CONCLUSIONS 

1.  The  DAFL  three-dimensional,  six-degrees-of-freedom  code  appears  to  be 
the  best  available  computer  model  for  predicting  projectile  penetration  in 
soil. 

2.  Many  items  in  the  DAFL  code  require  a more  precise  physical  description 
in  lieu  of  the  so-called  engineering  judgment  for  acceptable  understanding 
and  use  of  the  code.  Such  items  as  the  nose  shape  factor  and  shock  term 
are  examples  that  require  such  considerations. 

3.  When  simulating  AFATL  test  results,  with  the  code,  C^'s  varying  from  4 

to  2 at  high  velocity  for  dry  sand  gave  the  best  results.  The  inclusion  of 
a shock  term  for  initial  impact  would  most  likely  have  given  even  better 
predictive  results. 

4.  For  wet  sand,  C^  varying  from  six  at  zero  to  0.7  at  high  velocity  gave 

good  simulation.  Here  again,  an  initial  shock  term  would  have  enhanced  the 
predictive  results. 
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